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C — ' Abstract 

! ' These lectures are intended to give a pedagogical introduction to 

J> ■ the main current picture of the very early universe. After elementary 

reviews of general relativity and of the standard Big Bang model, the 
following subjects are discussed: inflation, the classical relativistic the- 
■ ory of cosmological perturbations and the generation of perturbations 

, from scalar field quantum fluctuations during inflation. 

o . 

*S : 1 Introduction 

The purpose of these lectures is to give an introduction to the present stan- 
dard picture of the early universe, which complements the older standard Big 
Bang model. These notes are intended for non-experts on this subject. They 
start with a very short introduction to General Relativity, on which modern 
cosmology is based, followed by an elementary review of the standard Big 
Bang model. We then discuss the limitations of this model and enter into 
the main subject of these lectures: inflation. 

Inflation was initially invented to solve some of the problems of the stan- 
dard Big Bang model and to get rid of unwanted relics generically predicted 
by high energy models. It turned out that inflation, as it was realized later, 
could also solve an additional puzzle of the standard model, that of the gener- 
ation of the cosmological perturbations. This welcome surprise put inflation 
on a rather firm footing, about twenty years ago. Twenty years later, infla- 
tion is still alive, in a stronger position than ever because its few competitors 
have been eliminated as new cosmological observations have accumulated 
during the last few years. 



1 



2 A few elements on general relativity and 
cosmology 

Modern cosmology is based on Einstein's theory of general relativity. It is 
thus useful, before discussing the early universe, to recall a few notions and 
useful formulas from this theory. Details can be found in standard textbooks 
on general relativity (see e.g. pQ). In the framework of general relativity, 
the spacetime geometry is defined by a metric, a symmetric tensor with two 
indices, whose components in a coordinate system (// = 0, 1, 2, 3) will be 
denoted g^ v . The square of the "distance" between two neighbouring points 
of spacetime is given by the expression 

ds 2 = g iiv dx il dx v . (1) 

We will use the signature (—,+,+,+). 

In a coordinate change — > x^, the new components of the metric are 
obtained by using the standard tensor transformation formulas, namely 

dx p dx a . . 

One can define a covariant derivative associated to this metric, denoted 
D^, whose action on a tensor with, for example, one covariant index and one 
contravariant index will be given by 

D X T\ = d x T» v + r» a T% - r%,T» ff (3) 

(a similar term must be added for each additional covariant or contravariant 
index), where the T are the Christoffel symbols (they are not tensors), defined 
by ^ 

= ^9 Xa {d^au + d v9iMJ - d a9ia ,) . (4) 

We have used the notation g^ v which corresponds, for the metric (and only 
for the metric), to the inverse of in a matricial sense, i.e. g^g™ = 5". 

The "curvature" of spacetime is characterized by the Riemann tensor, 
whose components can be expressed in terms of the Christoffel symbols ac- 
cording to the expression 

Rx,u p = 9,r{ u - ajv, + ti u t^ - v; u v p aX . (5) 

Einstein's equations relate the spacetime geometry to its matter content. 
The geometry appears in Einstein's equations via the Ricci tensor, defined 
by 

RpuV = R^av 1 (6) 
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and the scalar curvature, which is the trace of the Ricci tensor, i.e. 



R=<rR liV . (7) 

The matter enters Einstein's equations via the energy-momentum tensor, 
denoted T^, whose time/time component corresponds to the energy density, 
the time/space components to the momentum density and the space/space 
component to the stress tensor. Einstein's equations then read 

G^ u = Rp,, — -Rg^v = 8nGT^ u , (8) 

where the tensor G^ v is called the Einstein tensor. Since, by construction, 
the Einstein tensor satisfies the identity D^G 1 ^ = 0, any energy- momentum 
on the right-hand side of Einstein's equation must necessarily satisfy the 
relation 

DpT* = 0, (9) 

which can be interpreted as a generalization, in the context of a curved 
spacetime, of the familiar conservation laws for energy and momentum. 

The motion of a particule is described by its trajectory in spacetime, 
x M (A), where A is a parameter. A free particle, i.e. which does not feel any 
force (other than gravity), satisfies the geodesic equation, which reads 

fDX = 0, (10) 

where t M = dx^/dX is the vector field tangent to the trajectory (note that 
the geodesic equation written in this form assumes that the parameter A is 
afline). Equivalently, the geodesic can be rewritten as 

abf dx°_ 

d\ 2 pa d\ d\ 1 ' 

The geodesic equation applies both to 

• massive particles, in which case one usually takes as the parameter A 
the so-called proper time so that the corresponding tangent vector u p 
is normalized: g^ u u^u v — — 1; 

• massless particles, in particular the photon, in which case the tangent 
vector, usually denoted k p is light-like, i.e. g^k^ky = 0. 

Einstein's equations can also be obtained from a variational principle. 
The corresponding action reads 

S = ^£ / d^xsf^g (R - 2A) + J d A x^C mat . (12) 
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One can check that the variation of this action with respect to the metric 
g^ v , upon using the definition 

indeed gives Einstein's equations 

G lu , + Ag lu , = 8nGT IH ,. (14) 

This is a slight generalization of Einstein's equations (JBJ) that includes a 
cosmological constant A. It is worth noticing that the cosmo logical constant 
can also be interpreted as a particular energy-momentum tensor of the form 

2.1 Review of standard cosmology 

In this subsection, the foundations of modern cosmology are briefly recalled. 
They follow from Einstein's equations introduced above and from a few hy- 
potheses concerning spacetime and its matter content. One of the essential 
assumptions of cosmology (so far confirmed by observations) is to consider, as 
a first approximation, the universe as being homogeneous and isotropic. Note 
that these symmetries define implicitly a particular "slicing" of spacetime, 
the corresponding space-like hypersurfaces being homogeneous and isotropic. 
A different slicing of the same spacetime will give in general space-like hy- 
persurfaces that are not homogeneous and isotropic. 

The above hypothesis turns out to be very restrictive and the only metrics 
compatible with this requirement reduce to the so-called Robertson-Walker 
metrics, which read in an appropriate coordinate system 



ds 2 = -dt 2 + a 2 (t) 



llr " +r 2 (d9 2 + sin 2 Odcj> 2 ) 



1 — nr 



(15) 



with k = 0,-1,1 depending on the curvature of spatial hypersurfaces: re- 
spectively flat, elliptic or hyperbolic. 

The matter content compatible with the spacetime symmetries of homo- 
geneity and isotropy is necessarily described by an energy-momentum tensor 
of the form (in the same coordinate system as for the metric (|15p): 

T£ = Diag(-p(f),p(*),p(t),p(*))- (16) 

The quantity p corresponds to an energy density and P to a pressure. 

One can show that the so-called comoving particles, i.e. those particles 
whose spatial coordinates are constant in time, satisfy the geodesic equation 

(ED- 
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2.2 Friedmann-Lemaitre equations 



Substituting the Robertson- Walker metric (jl5j) in Einstein's equations 
one gets the so-called Friedmann-Lemaitre equations: 



-) = (17) 



d\ 2 8nGp K 

a J 3 a 

a 4ttG , 

(p + 3P). (1 



a 3 

An immediate consequence of these two equations is the continuity equation 

p + 3H(p + p)=0, (19) 

where H = a/a is the Hubble parameter. The continuity equation can be 
also obtained directly from the energy-momentum conservation D^T 1 ^ = 0, 
as mentioned before. 

In order to determine the cosmological evolution, it is easier to combine 
(|T7|) with (fT9"j) . Let us assume an equation of state for the cosmological 
matter of the form p = wp with w constant. This includes the two main 
types of matter that play an important role in cosmology: 

• gas of relativistic particles, w = 1/3; 

• non relativistic matter, w ~ 0. 

In these cases, the conservation equation (fT^j) can be integrated to give 

p oc a- 3(1+w) . (20) 
Substituting in (fTTjl . one finds, for k = 0, 

a 2 / a \ - 3 ( 1 + M ') 

3 gvrGpo - , (21) 

a A \a J 

where, by convention, the subscript '0' stands for present quantities. One 
thus finds a 2 oc a 2 ~ 3 ^ 1+w \ which gives for the evolution of the scale factor 

• in a universe dominated by non relativistic matter 

a(t) oc t 2/3 , (22) 

• and in a universe dominated by radiation 

a{t) oc t l/2 . (23) 
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One can also mention the case of a cosmological constant, which corre- 
sponds to an equation of state w = — 1 and thus implies an exponential 
evolution for the scale factor 

a{t) oc exp(Ht). (24) 

More generally, when several types of matter coexist with respectively 
P(i) = W(i)P(i), it is convenient to introduce the dimensionless parameters 

{t) ~ 3ff 2 ' ^ ' 

which express the present ratio of the energy density of some given species 
with respect to the so-called critical energy density p cr i t = 3Hq/ (8tcG), which 
corresponds to the total energy density for a flat universe. 
One can then rewrite the first Friedmann equation (fTTj) as 



/H\ 2 ^ /a\- 3 ( 1 +« , (.)) 



vi) =£M^ +n -U) • < 26 > 

with Q K = —k/o^Hq, which implies that the cosmological parameters must 
satisfy the consistency relation 

£%)+0 K = l. (27) 

i 

As for the second Friedmann equation (fT%|) . it implies 

= -^^ )(1+Ww) - (28) 
Present cosmological observations yield for the various parameters 

• Baryons: Qb — 0.05, 

• Dark matter: ~ 0.25, 

• Dark energy (compatible with a cosmological constant): — 0.7, 

• Photons: ^ 7 ~5x 10 -5 . 

Observations have not detected so far any deviation from flatness. Radiation 
is very subdominant today but extrapolating backwards in time, radiation 
was dominant in the past since its energy density scales as p 7 oc a -4 in 
contrast with non relativistic matter (p m oc a -3 ). Moreover, since the present 
matter content seems dominated by dark energy similar to a cosmological 
constant (w\ = —1), this indicates that our universe is presently accelerating. 
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2.3 The cosmological redshift 

An important consequence of the expansion of the universe is the cosmological 
redshift of photons. This is in fact how the expansion of the universe was 
discovered initially. 

Let us consider two light signals emitted by a comoving object at two 
successive instants t e and t e + St e , then received later at respectively t Q and 
t + 8t a by a (comoving) observer. One can always set the observer at the 
center of the coordinate system. All light trajectories reaching the observer 
are then radial and one can write, using (fTHjl 

dr f*° dt 



o y/l — nr 2 Jt e a(t) 

The left-hand side being identical for the two successive trajectories, the 
right-hand side must vanish, which yields 

5t„ 8t P 

— = 0. (30) 

a Q a e 

This implies for the frequencies measured at emission and at reception a 
redshift given by 

l + z = ^l = ^. (3i) 



2.4 Thermal history of the universe 

To go beyond a simply geometrical description of cosmology, it is very fruitful 
to apply thermodynamics to the matter content of the universe. One can 
then define a temperature T for the cosmological photons, not only when 
they are strongly interacting with ordinary matter but also after they have 
decoupled because, with the expansion, the thermal distribution for the gas 
of photons is unchanged except for a global rescaling of the temperature so 
that T essentially evolves as 

T(t) oc (32) 

This means that, going backwards in time, the universe was hotter and hotter. 
This is the essence of the hot Big Bang scenario. 

As the universe evolves, the reaction rates between the various species 
are modified. A detailed analysis of these changes enables to reconstruct 
the past thermal history of the universe. Two events in particular play an 
essential role because of their observational consequences: 
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• Primordial nucleosynthesis 

Nucleosynthesis occured at a temperature around 0.1 MeV, when the 
average kinetic energy became sufficiently low so that nuclear bind- 
ing was possible. Protons and neutrons could then combine, which 
lead to the production of light elements, such that Helium, Deuterium, 
Lithium, etc... Within the observational uncertainties, this scenario is 
remarkably confirmed by the present measurements. 

• Decoupling of baryons and photons (or last scattering) 

A more recent event is the so-called "recombination" of nuclei and 
electrons to form atoms. This occured at a temperature of the order 
of the eV. Free electrons thus almost disappeared, which entailed an 
effective decoupling of the cosmological photons and ordinary matter. 
What we see today as the Cosmic Microwave Background (CMB) is 
made of the fossil photons, which interacted for the last time with 
matter at the last scattering epoch. The CMB represents a remarkable 
observational tool for analysing the perturbations of the early universe, 
as well as for measuring the cosmological parameters introduced above. 

2.5 Puzzles of the standard Big Bang model 

The standard Big Bang model has encountered remarkable successes, in par- 
ticular with the nucleosynthesis scenario and the prediction of the CMB, and 
it remains today a cornerstone in our understanding of the present and past 
universe. However, a few intriguing facts remain unexplained in the strict 
scenario of the standard Big Bang model and seem to necessitate a larger 
framework. We review below the main problems: 

• Homogeneity problem 

A first question is why the approximation of homogeneity and isotropy 
turns out to be so good. Indeed, inhomogeneities are unstable, because 
of gravitation, and they tend to grow with time. It can be verified 
for instance with the CMB that inhomogeneities were much smaller 
at the last scattering epoch than today. One thus expects that these 
homogeneities were still smaller further back in time. How to explain 
a universe so smooth in its past ? 

• Flatness problem 
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Another puzzle lies in the (spatial) flatness of our universe. Indeed, 
Friedmann's equation (JT7j) implies 

In standard cosmology, the scale factor behaves like a ~ t p with p < 1 
(p = 1/2 for radiation and p = 2/3 for non-relativistic matter). As a 
consequence, (aH)~ 2 grows with time and \Q — 1| must thus diverge 
with time. Therefore, in the context of the standard model, the quasi- 
flatness observed today requires an extreme fine-tuning of Q near 1 in 
the early universe. 

• Horizon problem 

One of the most fundamental problems in standard cosmology is cer- 
tainly the horizon problem. The (particle) horizon is the maximal dis- 
tance that can be covered by a light ray. For a light-like radial trajec- 
tory dr = a(t)dt and the horizon is thus given by 

d ^ = a ^l',W) =a{t)t ^T 1 ' (34) 

where the last equality is obtained by assuming a(t) ~ t q and tj is some 
initial time. 

In standard cosmology (q < 1), the integral converges in the limit = 
and the horizon has a finite size, of the order of the so-called Hubble 
radius H~ x : 

d H (t) = -r—H- 1 . (35) 
1 - q 

It also useful to consider the comoving Hubble radius, (aif) -1 , which 
represents the fraction of comoving space in causal contact. One finds 
that it grows with time, which means that the fraction of the universe in 
causal contact increases with time in the context of standard cosmology. 
But the CMB tells us that the Universe was quasi-homogeneous at the 
time of last scattering on a scale encompassing many regions a priori 
causally independent. How to explain this ? 

A solution to the horizon problem and to the other puzzles is provided 
by the inflationary scenario, which we will examine in the next section. The 
basic idea is to invert the behaviour of the comoving Hubble radius, that is to 
make him decrease sufficiently in the very early universe. The corresponding 
condition is that 

a > 0, (36) 
i.e. that the Universe must undergo a phase of acceleration. 
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Figure 1: Evolution of the comoving Hubble radius Xh = (aif) -1 : during 
standard cosmology, Xh increases (continuous lines), whereas during inflation 
\h shrinks (dashed lines). 
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3 Inflation 



The broadest definition of inflation is that it corresponds to a phase of ac- 
celeration of the universe, 

a > 0. (37) 

In this broad sense, the current cosmological observations, if correctly inter- 
preted, mean that our present universe is undergoing an inflationary phase. 
We are however interested here in an inflationary phase taking place in the 
very early universe, with different energy scales. 

The Friedmann equations (JT7j) tell us that one can get acceleration only 
if the equation of state satisfies the condition 

P < ~\p, (38) 

condition which looks at first view rather exotic. 

A very simple example giving such an equation of state is a cosmological 
constant, corresponding to a cosmological fluid with the equation of state 

P = -p. (39) 

However, a strict cosmological constant leads to exponential inflation forever 
which cannot be followed by a radiation or matter era. Another possibility 
is a scalar field, which we discuss now in some details. 

3.1 Cosmological scalar fields 

The dynamics of a scalar field coupled to gravity is governed by the action 
SV = y d 4 x^9 [~\d^d^ - V(0)) . (40) 

The corresponding energy-momentum tensor, which can be derived using 
(|T3|). is given by 

= <VR0 - 9ia , f~d CT <R0 + V(<j>j) . (41) 

If one assumes the geometry, and thus the matter, to be homogeneous and 
isotropic, then the energy-momentum tensor reduces to the perfect fluid form 
with the energy density 

P = -T° = ^ + V(<J ) ) } (42) 
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where one recognizes the sum of a kinetic energy and of a potential energy, 
and the pressure 

P=\^-V(<f>). (43) 

The equation of motion for the scalar field is the Klein-Gordon equation, 
which is obtained by taking the variation of the above action (140)1 with respect 
to the scalar field and which reads 

D^Drf = V, (44) 

in general and 

<j> + 3H<j) + V = (45) 

in the particular case of a FLRW (Friedmann-Lemaitre-Robertson- Walker) 
universe. 

The system of equations governing the dynamics of the scalar field and 
of the geometry in a FLRW universe is thus given by 

# 2 = ^(^ 2 + W), (46) 

4> + 3H0 + V' = 0, (47) 
H = -4ttG0 2 . (48) 

The last equation can be derived from the first two and is therefore redun- 
dant. 



3.2 The slow-roll regime 



The dynamical system (j46M48j) does not always give an accelerated expansion 
but it does so in the so-called slow-roll regime when the potential energy of 
the scalar field dominates over its kinetic energy. 

More specifically, the so-called slow roll approximation consists in neglect- 
ing the kinetic energy of the scalar field , 2 /2 in (jlBj) and the acceleration 
in the Klein- Gordon equation (|47|) . One then gets the simplified system 

H> ^ "-fv, (49) 

2,H4> + V ~ 0. (50) 

Let us now examine in which regime this approximation is valid. From (|5()j) . 
the velocity of the scalar field is given by 

*"-w (51) 
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Substituting this relation in the condition <f) 2 /2 <C V yields the requirement: 



where we have introduced the reduced Planck mass 

m P = (53) 

Similarly, the time derivative of ([51)1. using the time derivative of (@HJ), gives, 
after substitution in <fi <C V, the condition 

V" 

T) V = m P — < 1. (54) 

In summary, the slow-roll approximation is valid when the two conditions 
€y, Vv 1 are satisfied, which means that the slope and the curvature of 
the potential, in Planck units, must be sufficiently small. 

3.3 Number of e-folds 

When working with a specific inflationary model, it is important to be able 
to relate the cosmological scales observed at the present time with the scales 
during inflation. For this purpose, one usally introduces the number of e- 
foldings before the end of inflation, denoted N, and simply defined by 

iV = ln^, (55) 
a 

where a en d is the value of the scale factor at the end of inflation and a is a 
fiducial value for the scale factor during inflation. By definition, N decreases 
during the inflationary phase and reaches zero at its end. In the slow-roll 
approximation, it is possible to express N as a function of the scalar field. 
Since dN = — dlna = —Hdt = —(H/<j))d<j), one easily finds, using (|51|) and 
gUl), that 

*W=7, (56) 

Given an explicit potential V(4>), one can in principle integrate the above 
expression to obtain N in terms of 0. This will be illustrated below for some 
specific models. 

Let us now discuss the link between iV and the present cosmological scales. 
Let us consider a given scale characterized by its comoving wavenumber 
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k = 2n/\. This scale crossed outside the Hubble radius, during inflation, at 
an instant t*(k) defined by 

k = a(QH(Q. (57) 

To get a rough estimate of the number of e-foldings of inflation that are 
needed to solve the horizon problem, let us first ignore the transition from 
a radiation era to a matter era and assume for simplicity that the inflation- 
ary phase was followed instantaneously by a radiation phase that has lasted 
until now. During the radiation phase, the comoving Hubble radius (aH)^ 1 
increases like a. In order to solve the horizon problem, the increase of the co- 
moving Hubble radius during the standard evolution must be compensated 
by at least a decrease of the same amount during inflation. Since the co- 
moving Hubble radius roughly scales like a -1 during inflation, the minimum 
amount of inflation is simply given by the number of e-folds between the end 
of inflation and today Ha /a end ) = \n(T end /T ) ~ ln(10 29 (T end /10 16 GeV)), 
i.e. around 60 e-folds for a temperature T ~ 10 16 Gev at the beginning of the 
radiation era. As we will see later, this energy scale is typical of inflation in 
the simplest models. 

This determines roughly the number of e-folds N(k ) between the mo- 
ment when the scale corresponding to our present Hubble radius k = a H 
exited the Hubble radius during inflation and the end of inflation. The other 
lengthscales of cosmological interest are smaller than k^ 1 and therefore exited 
the Hubble radius during inflation after the scale ko, whereas they entered 
the Hubble radius during the standard cosmological phase (either in the ra- 
diation era for the smaller scales or in the matter era for the larger scales) 
before the scale k . 

A more detailed calculation, which distinguishes between the energy scales 
at the end of inflation and after the reheating, gives for the number of e-folds 
between the exit of the mode k and the end of inflation 

7 y-1/4 yl/4 , 1/4 

= 62 " ' n ^ + WkeV + k Wi + 3 k Wf <58) 

v end K end 

Since the smallest scale of cosmological relevance is of the order of 1 Mpc, 
the range of cosmological scales covers about 9 e-folds. 

The above number of e-folds is altered if one changes the thermal his- 
tory of the universe between inflation and the present time by including for 
instance a period of so-called thermal inflation. 
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3.4 Power- law potentials 



It is now time to illustrate all the points discussed above with some specific 
potential. We consider first the case of power-law monomial potentials, of 
the form 



v{4>) 



A 

—5 

V 



-nil 



nip j 



(59) 



which have been abundantly used in the literature. In particular, the above 
potential includes the case of a free massive scalar field, V(4>) = m 2 0/2. The 
slow-roll conditions e < 1 and rj <C 1 both imply 



~> p nip, 



(60) 



which means that the scalar field amplitude must be above the Planck mass 
during inflation. 

After subsituting the potential (|5T)j) into the slow-roll equations of motion 
one can integrate them explicitly to get 



,2-« 



for p ^ 4 and 



exp 



4-p 
4A 



— nip (t-ti) 



nip{t - U) 



(61) 



(62) 



for p — 4. 

One can also express the scale factor as a function of the scalar field (and 
thus as a function of time by substituting the above expression for </>(£)) by 
using d\na/d(f) = H/<f> ~ —<f)/(pmp). One finds 



a = a end exp 



2pmp 



(63) 



Defining the end of inflation by ey = 1, one gets 4> enc i = pm P /\/2 and the 
number of e-folds is thus given by 



N{<f>) 



2pm 



P 



2 4 - 



This can be inverted, so that 



4>{N) ~ ^2Npm P , 



(64) 



(65) 



where we have ignored the second term of the right hand side of (I64j) . in 
agreement with the condition (|6()jl. 
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3.5 Exponential potential 



If one considers a potential of the form 

V = V exp (-* 



12 
q m P 



(66) 



then it is possible to find an exact solution (i.e. valid beyond the slow-roll 
approximation) of the system (|46H48j) . with a power-law scale factor, i.e. 

(67) 



a(t) oc t q . 

The evolution of the scalar field is given by the expression 

t 



4>(t) 



2q rap In 



Vn 



q(3q — 1) mp 



(6* 



Note that one recovers the slow-roll approximation in the limit g > 1, since 
the slow-roll parameters are given by ey = 1/q and r/v = 2/q. 



3.6 Brief history of the inflationary models 

Let us now try to summarize in a few lines the history of inflationary mod- 
els. The first model of inflation is usually traced back to Alan Guth [2] in 
1981, although one can see a precursor in the model of Alexei Starobinsky [3]. 
Guth's model, which is named today old inflation is based on a first-order 
phase transition, from a false vacuum with non zero energy, which generates 
an exponential inflationary phase, into a true vacuum with zero energy den- 
sity. The true vacuum phase appears in the shape of bubbles via quantum 
tunneling. The problem with this inflationary model is that, in order to get 
sufficient inflation to solve the problems of the standard model mentioned 
earlier, the nucleation rate must be sufficiently small; but, then, the bubbles 
never coalesce because the space that separates the bubbles undergoes infla- 
tion and expands too rapidly. Therefore, the first model of inflation is not 
phenomenologically viable. 

After this first and unsuccessful attempt, a new generation of inflationary 
models appeared, usually denoted new inflation models jlj. They rely on a 
second order phase transition, based on thermal corrections of the effective 
potential and thus assume that the scalar field is in thermal equilibrium. 

This hypothesis of thermal equilibrium was given up in the third genera- 
tion of models, initiated by Andrei Linde, and whose generic name is chaotic 
inflation [5J. This allows to use extremely simple potentials, quadratic or 
quartic, which lead to inflationary phases when the scalar field is displaced 
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from the origin with values of the order of several Planck masses. This is 
sometimes considered to be problematic from a particle physics point of view, 
as discussed briefly later. 

During the last few years, there has been a revival of the inflation model 
building based on high energy theories, in particular in the context of super- 
symmetry. In these models, the value of the scalar field is much smaller than 
the Planck mass. 

3.7 The inflationary "zoology" 

There exist many models of inflation. As far as single-field models are con- 
cerned (or at least effectively single field during inflation, the hybrid models 
requiring a second field to end inflation as discussed below), it is convenient 
to regroup them into three broad categories: 

• Large field models (0 < rj < e) 

The scalar field is displaced from its stable minimum by A0 ~ mp. 
This includes the so-called chaotic type models with monomial poten- 
tials 

\/(0) = A 4 ^ P , (69) 

or the exponential potential 

V(0) = A 4 exp (0//i) , (70) 
which have already been discussed. 

This category of models is widely used in the literature because of the 
computational simplicity. But they are not always considered to be 
models which can be well motivated by particle physics. The reason is 
the following. The generic potential for a scalar field will contain an 
infinite number of terms, 

V{4>) = V + V0 2 + f 3 + f 4 + E <V. (71) 

where the non-renormalizable (d > 4) couplings are a priori of order 
1. When the scalar field is of order of a few Planck masses, one has 
no control on the form of the potential, and all the non-normalizable 
terms must be taken into account in principle. 

In order to work with more specific forms for the potential, inflationary 
model builders tend to concentrate on models where the scalar field 
amplitude is small with respect to the Planck mass, as those discussed 
just below. 
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Figure 2: Schematic potential for t he three main categories of inflationary 
models: a. chaotic models b. symmetry breaking models; c. hybrid models 

• Small field models (77 < < e) 

In this type of models, the scalar field is rolling away from an unstable 
maximum of the potential. This is a characteristic feature of sponta- 
neous symmetry breaking. A typical potential is 



which can be interpreted as the lowest-order term in a Taylor expansion 
about the origin. Historically, this potential shape appeared in the so- 
called 'new inflation' scenario. 

A particular feature of these models is that tensor modes are much 
more suppressed with respect to scalar modes than in the large-field 
models, as it will be shown later. 

• Hybrid models (0 < e < 77) 

This category of models, which appeared rather recently, relies on the 
presence of two scalar fields: one plays the traditional role of the infla- 
ton, while the other is necessary to end inflation. 

As an illustration, let us consider the original model of hybrid inflation 
based on the potential 



For values of the field (f> larger than the critical value <p c = XM 2 /X', 
the potential for ip has its minimum at ip — 0. This is the case during 
inflation, ip is thus trapped in this minimum ip = 0, so that the effective 
potential for the scalar field 0, which plays the role of the inflaton, is 
given by 




(72) 




(73) 



^//(0) = V o + -m 2 2 , 



(74) 
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Figure 3: Typical potential V((j),ip) for hybrid inflation. 
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with Vq = AM 4 /4. During the inflationary phase, the field <fi slow-rolls 
until it reaches the critical value (f> c . The shape of the potential for if) 
is then modified and new minima appear in if> = ±M. if) will thus roll 
down into one of these new minima and, as a consequence, inflation 
will end. 

As far as inflation is concerned, hybrid inflation scenarios correspond ef- 
fectively to single-field models with a potential characterized by V"((f>) > 
and < e < r\. A typical potential is 



Once more, this potential can be seen as the lowest order in a Taylor 
expansion about the origin. 

In the case of hybrid models, the value 4>n of the scalar field as a func- 
tion of the number of e-folds before the end of inflation is not deter- 
mined by the above potential and, therefore, ((f) n/ A 4 ) can be considered 
freely adjustable parameter. 

Many more details on inflationary models can be found in e.g. EJ El- 

4 The theory of cosmological perturbations 

So far, we have concentrated our attention on strictly homogeneous and 
isotropic aspects of cosmology. Of course, this idealized version, although 
extremely useful, is not sufficient to account for real cosmology and it is now 
time to turn to the study of deviations from homogeneity and isotropy. 

In cosmology, inhomogeneities grow because of the attractive nature of 
gravity, which implies that inhomogeneities were much smaller in the past. 
As a consequence, for most of their evolution, inhomogeneities can be treated 
as linear perturbations. The linear treatment ceases to be valid on small scales 
in our recent past, hence the difficulty to reconstruct the primordial inho- 
mogeneities from large-scale structure, but it is quite adequate to describe 
the fluctuations of the CMB at the time of last scattering. This is the rea- 
son why the CMB is currently the best observational probe of primordial 
inhomogeneities. For more details on the relativistic theory of cosmological 
perturbations, which will be briefly introduced in this chapter, the reader is 
invited to read the standard reviews jTU] in the literature. 

From now on, we will be mostly working with the conformal time rj, 
instead of the cosmic time t. The conformal time is defined as 




(75) 




(76) 
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so that the (spatially flat) FLRW metric takes the remarkably simple form 



ds 2 = a 2 (r]) -drj 2 + 5 ij dx i dx j . (77) 
4.1 Perturbations of the geometry 

Let us start with the linear perturbations of the geometry. The most general 
linear perturbation of the FLRW metric can be expressed as 

ds 2 = a 2 {-(1 + 2A)drf + 2B i dx i dr) + (<% + h^) dx*dx j } , (78) 

where we have considered only the spatially flat FLRW metric. 

We have introduced a time plus space decomposition of the perturbations. 
The indices i, j stand for spatial indices and the perturbed quantities defined 
in (J78j) can be seen as three-dimensional tensors, for which the indices can 
be lowered (or raised) by the spatial metric <5y (or its inverse). 

It is very convenient to separate the perturbations into three categories, 
the so-called "scalar" , "vector" and "tensor" modes. For example, a spatial 
vector field B l can be decomposed uniquely into a longitudinal part and a 
transverse part, 

B % = ViB + B i} ViB 1 = 0, (79) 

where the longitudinal part is curl-free and can thus be expressed as a gra- 
dient, and the transverse part is divergenceless. One thus gets one "scalar" 
mode, B , and two "vector" modes B l (the index i takes three values but the 
divergenceless condition implies that only two components are independent. 

A similar procedure applies to the symmetric tensor hij, which can be 
decomposed as 

Kj = 2C<% + 2ViVjE + 2V {i E j) + E ih (80) 

with E %1 transverse and traceless (TT), i.e. VjE 13 = (transverse) and 
E t] 5ij = (traceless), and Ei transverse. The parentheses around the indices 
denote symmetrization, namely 2V(j-EV) = Vj-£y + V ' jE^. We have thus 
defined two scalar modes, C and E, two vector modes, Ei, and two tensor 
modes, Eij. 

In the following, we will be mainly interested in the metric with only 
scalar perturbations, since scalar modes are the most relevant for cosmology 
and they can be treated independently of the vector and tensor modes. In 
matrix notation, the perturbed metric will thus be of the form 



g/iu = a 2 



1 + 2A) V Z B 
VjB {(1 + 2C)5ij + 2ViVjE} 
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After the description of the perturbed geometry, we turn to the perturbations 
of the matter in the next subsection. 



4.2 Perturbations of the matter 



Quite generally, the perturbed energy-momentum tensor can be written in 
the form 

~{p + Sp) qj 



Tff 



-q* + (p + P)B i (P + 8P)5) + ir) 



(82) 



where is the momentum density and ttj is the anisotropic stress tensor, 
which is traceless, i.e. = 0. One can then decompose these tensors into 
scalar, vector and tensor components, as before for the metric components, 
so that 

qi = V iq + q h V^ = 0, (83) 

and 



7T 



V i V i 7T--5 ii V fc V ; 



7T- 



-2V 



V fc 7T fc = 0, 



0. 



7T 



(84) 



4.2.1 Fluid 

A widely used description for matter in cosmology is that of a fluid. Its 
homogeneous part is described by the energy-momentum tensor of a perfect 
fluid, as seen earlier, while its perturbed part can be expressed as 

5T °= -Sp, (85) 
6T*= -(p + p)v\ 5Tf = {p + p){v l + B l ) (86) 
5T;= SP5) + ni, (87) 

with 7r^ = as before and where v l is the three-dimensional fluid velocity 
defined by 

6u { = -v\ (88) 
a 

It is also possible to separate this perturbed energy-momentum tensor into 
scalar, vector and tensor parts by using decompositions for and 7ry similar 
to (JHl and (JHU). 
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4.2.2 Scalar field 



Another type of matter, especially useful in the context of inflation, is a scalar 
field. The homogeneous description has already been given earlier and the 
perturbed expression for the energy-momentum tensor follows immediately 
from (j41j) . taking into account the metric perturbations as well. One finds 

-a 2 5p = -<f>'5<// - a 2 V'5(f) + <j)' 2 A, (89) 
a 2 qi = -0'di64>, (90) 
-5) (a 2 V'5(t) + <p' 2 A - <j>'5<jt) . (91) 

The last equation shows that, for a scalar field, there is no anisotropic stress 
in the energy-momentum tensor. 

4.3 Gauge transformations 

It is worth noticing that there is a fundamental difference between the per- 
turbations in general relativity and the perturbations in other field theories 
where the underlying spacetime is fixed. In the latter case, one can define 
the perturbation of a given field </> as 

5<P(p)=<P(p)-0(p), (92) 

where 4> is the unperturbed field and p is any point of the spacetime. 

In the context of general relativity, spacetime is no longer a frozen back- 
ground but must also be perturbed if matter is perturbed. As a consequence, 
the above definition does not make sense since the perturbed quantity lives 
in the perturbed spacetime Ai, whereas the unperturbed quantity lives in 
another spacetime: the unperturbed spacetime of reference, which we will 
denote Ai. In order to use a definition similar to ([92)1 . one must introduce a 
one-to-one identification, i, between the points of Ai and those of Ai. The 
perturbation of the field can then be defined as 

5<P(p) = MP))-m, (93) 

where p is a point of Ai. 

However, the identification i is not uniquely defined, and therefore the 
definition of the perturbation depends on the particular choice for i: two 
different identifications, i\ and ii say, lead to two different definitions for the 
perturbations. This ambiguity can be related to the freedom of choice of the 
coordinate system. Indeed, if one is given a coordinate system in At, one can 
transport it into At via the identification. L\ and 12 thus define two different 



a 2 5T ° = 
a 2 5T° = 
a 2 bT] = 
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coordinate systems in A4, and in this respect, a change of identification can 
be seen as a change of coordinates in M.. 

The metric perturbations, introduced in (|7Hj) . are modified in a coordinate 
transformation of the form 

x a ^x a +e, r = (£°,e). (94) 

It can be shown that the change of the metric components can be expressed 
as 

A{6g lw ) = -2D Qt ^ ) . (95) 

where A represents the variation, due the coordinate transformation, at the 
same old and new coordinates (and thus at different physical points). The 
above variation can be decomposed into individual variations for the various 
components of the metric defined earlier. One finds 

AA = -t°'-TiZ (96) 
AB t = V^-i'i (97) 
Ahij = -2 (V (i £i) - ■ (98) 

The effect of a coordinate transformation can also be decomposed along 
the scalar, vector and tensor sectors introduced earlier. The generator of 
the coordinate transformation can be written as 

r = (£°,v 4 e+r), (99) 

with transverse. This shows explicitly that £ a contains two scalar compo- 
nents, £° and £, and two vector components, £ . The transformations (|96H98|) 
are then decomposed into : 

a -> A-f-n^ 
c -> c-nz° 

E ^ E-£ (100) 

if -> if-f' 
^ -» 

The tensor perturbations remain unchanged since £ a does not contain any 
tensor component. The matter perturbations, either in the fluid description 
or in the scalar field description, follow similar transformation laws in a 
coordinate change. 



24 



In order to study the physically relevant modes and not spurious modes 
due to coordinate ambiguities, two strategies can a priori be envisaged. The 
first consists in working from the start in a specific gauge. A familiar choice 
in the literature on cosmological perturbations is the longitudinal gauge (also 
called conformal Newton gauge), which imposes 

B L = 0, E L = 0. (101) 

The second approach consists in defining gauge-invariant variables, i.e. 
variables that are left unchanged under a coordinate transformation. For the 
scalar metric perturbations, we start with four quantities (A, B, C and E) 
and we can use two gauge transformations (£° and £). This implies that the 
scalar metric perturbations must be described by two independent gauge- 
invariant quantities. Two such quantities are 

$ = A + (B - E')' + H(B - E') (102) 

and 

* = -C -H{B - E'), (103) 

as it can be checked by considering the explicit transformations in (jlOOj) . It 
turns out that, in the longitudinal gauge, the remaining scalar perturbations 
A L and Cl are numerically equivalent to the gauge-invariant quantities just 
defined <3> and — 

In practice, one can combine the two strategies by doing explicit calcu- 
lations in a given gauge and then by relating the quantities defined in this 
gauge to some gauge-invariant variables. It is then possible to translate the 
results in any other gauge. In the rest of these lectures, we will use the 
longitudinal gauge. 



4.4 The perturbed Einstein equations 

After having defined the metric and the matter perturbations, we can now 
relate them via the perturbed Einstein equations. We will consider here 
explicitly only the scalar sector, which is the most complicated but also the 
most interesting for cosmological applications. 

Starting from the perturbed metric (|81)L one can compute the components 
of Einstein's tensor at linear order. In the longitudinal gauge, i.e. with 
Bl = El = 0, one finds 



l a 2 



r 3H 2 A L -3HCl + V 2 C L ] (104) 



(<5G*°) l = ^Vi[-WA L + C[] (105) 
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( 5G l)r = \^j(-CL-A L ) + \\-2C^-AHCi + V 2 C L 

+27M' L + V 2 A L + 2 (2W + H 2 ) A L ] 8). (106) 



Combining with the perturbations of the energy-momentum tensor given 
in (|82jl. the perturbed Einstein equations yield, in the longitudinal gauge, the 
following relations: the energy constraint (from (jl()4j) ) 

3H 2 A L - 3HCL + V 2 C L = -47rGa 2 5p L , (107) 

the momentum constraint (from ([105)) ) 

C[ - HA L = AttG a 2 q L} (108) 
the "anisotropy constraint" (from the traceless part of ([106)) ) 

- A L - C L = 8ttG a 2 n L , (109) 

and finally 

C'l + 2HC' h - HA' L - (2W + H 2 )A L - ^V 2 (A L + C L ) = -4nG5P L , (110) 

obtained from the trace of ([106)1 . 

The combination of the energy and momentum constraints gives the useful 
relation 

V 2 ^ = AnGa 2 (5p L - 3Hq L ) = AnGa 2 5p c , (111) 

where we have introduced the comoving energy density perturbation Sp c : 
this gauge-invariant quantity corresponds, according to its definition, to the 
energy density perturbation measured in comoving gauges characterized by 
5Tj° = qi = 0. We have also replaced Cl by — \&. Note that the above equation 
is quite similar to the Newtonian Poisson equation, but with quantities whose 
natural interpretation is given in different gauges. 



4.5 Equations for the matter 

As mentioned earlier, a consequence of Einstein's equations is that the to- 
tal energy-momentum tensor is covariantly conserved (see Eq. (j^J). For a 
fluid, the conservation of the energy-momentum tensor leads to a continuity 
equation that generalizes the continuity equation of fluid mechanics, and a 
momentum conservation equation that generalizes the Euler equation. In the 
case of a single fluid, combinations of the perturbed Einstein equations ob- 
tained in the previous subsection lead necessarily to the perturbed continuity 
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and Euler equations for the fluid. In the case of several non-interacting flu- 
ids, however, one must impose separately the covariant conservation of each 
energy-momentum tensor: this is not a consequence of Einstein's equations, 
which impose only the conservation of the total energy-momentum tensor. 

In the perspective to deal with several cosmological fluids, it is therefore 
useful to write the perturbation equations, satisfied by a given fluid, that 
follow only from the conservation of the corresponding energy-momentum 
tensor, and independently of Einstein's equations. 

The continuity equation can be obtained by perturbing u^D v T^ = 0. One 
finds, in any gauge, 

bp' + m (bp + 5P) + (p + P) (SC + V 2 E' + V 2 v) =0. (112) 

Dividing by p, this can be reexpressed in terms of the density contrast 5 = 
5p/p: 

5' + m \ 5 -f- - w ) 5 + (1 + w) (v 2 v + 3C + V 2 E') = 0, (113) 



bp 

where w = p/p (w is not necessarily constant here). The perturbed Euler 
equation is derived from the spatial projection of 5(D u Tp = 0. This gives 

(v + B)' + (l - 3c 2 ) H(v + B) + A+ + 3(p ^ p) V 2 7T = 0, (114) 

where c s is the sound speed, which is related to the time derivatives of the 
background energy density and pressure: 

4 - y ("5) 

There are as many systems of equations (jll3H114|) as the number of fluids. If 
the fluids are interacting, one must add an interaction term on the right-hand 
side of the Euler equations. 

Finally, let us stress that the fluid description is not always an adequate 
approximation for cosmological matter. A typical example is the photons 
during and after recombination: their mean free path becomes so large that 
they must be treated as a gas, which requires the use of the Boltzmann 
equation (see e.g. 11 for a presentation of the Boltzmann equation in the 
cosmological context). 
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4.6 Initial conditions for standard cosmology 

The notion of initial conditions depends in general on the context, since 
the initial conditions for a given period in the history of the universe can 
be seen as the final conditions of the previous phase. In cosmology, "initial 
conditions" usually refer to the state of the perturbations during the radiation 
dominated era (of standard cosmology) and on wavelengths larger than the 
Hubble radius. 

Let us first address in details the question of initial conditions in the 
simple case of a single perfect fluid, radiation, with equation of state p = p/3 
(which gives c\ — w — 1/3). The four key equations are the continuity, Euler, 
Poisson and anisotropy equations, respectively Eqs ()113|) . ()114|) . (jlll|) and 
()109|) . In terms of the Fourier components, 

Q(k) = J (^ e ~^(^ ( 116 ) 
and of the dimensionless quantity 

x = ki] (117) 

(during the radiation dominated era 7i = 1/r]), the four equations can be 
rewritten as 

^-iv + 4^ = 0, (118) 
dx 6 dx 

% + \* + A = 0, (119) 
x 2 C=\(s--V\ (120) 



2 V x 

C = -A. (121) 

We have introduced the quantity V = kv, which has the dimension of a 
velocity. Since we are interested in perturbations with wavelength larger 
than the Hubble radius, i.e. such that x = k/TC <C 1, it is useful to consider 
a Taylor expansion for the various perturbations, for instance 

V = V (0) + xV (1) + y V (2) + . . . (122) 

One then substitutes these Taylor expansions into the above system of equa- 
tions. In particular, the Poisson equation ()120|) gives 



V (0) = 0, (123) 
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in order to avoid a divergence, as well as 



V). 



(124) 



4 



The Euler equation (jll9|) then gives 



A (o) 



(125) 
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The conclusion is that the initial conditions for each Fourier mode are deter- 
mined by a single quantity, e.g. S^°\ the other quantities being related via 
the above constraints. 

In general, one must consider several cosmological fluids. Typically, the 
"initial" or "primordial" perturbations are defined deep in the radiation era 
but at temperatures low enough, i.e. after nucleosynthesis, so that the main 
cosmological components reduce to the usual photons, baryons, neutrinos and 
cold dark matter (CDM). The system (jll8H12l|) must thus be generalized to 
include a continuity equation and a Euler equation for each fluid. The above 
various cosmological species can be characterized by their number density, 
n x , and their energy density px- In a multi-fluid system, it is useful to 
distinguish adiabatic and isocurvature perturbations. 

The adiabatic mode is defined as a perturbation affecting all the cosmo- 
logical species such that the relative ratios in the number densities remain 
unperturbed, i.e. such that 



It is associated with a curvature perturbation, via Einstein's equations, since 
there is a global perturbation of the matter content. This is why the adiabatic 
perturbation is also called curvature perturbation. In terms of the energy 
density contrasts, the adiabatic perturbation is characterized by the relations 



They follow directly from the prescription ()126|) . each coefficient depending 
on the equation of state of the particuler species. 

Since there are several cosmological species, it is also possible to perturb 
the matter components without perturbing the geometry. This corresponds 
to isocurvature perturbations, characterized by variations in the particle num- 
ber ratios but with vanishing curvature perturbation. The variation in the 
relative particle number densities between two species can be quantified by 
the so-called entropy perturbation 



5 {n x /n Y ) = 0. 



(126) 




(127) 
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When the equation of state for a given species is such that w = p/p = Const, 
then one can reexpress the entropy perturbation in terms of the density 
contrast, in the form 

sa* - T~~f~ - tt~- (129) 

1 + w A l + w B 

It is convenient to choose a species of reference, for instance the photons, 
and to define the entropy perturbations of the other species relative to it: 

(130) 
(131) 
(132) 

thus define respectively the baryon isocurvature mode, the CDM isocurvature 
mode, the neutrino isocurvature mode. In terms of the entropy perturbations, 
the adiabatic mode is obviously characterized by Sb = S c = S u = 0. 

In summary, we can decompose a general perturbation, described by four 
density contrasts, into one adiabatic mode and three isocurvature mode. In 
fact, the problem is slightly more complicated because the evolution of the 
initial velocity fields. For a single fluid, we have seen that the velocity field 
is not an independent initial condition but depends on the density contrast 
so that there is no divergence backwards in time. In the case of the four 
species mentioned above, there remains however one arbitrary relative veloc- 
ity between the species, which gives an additional mode, usually named the 
neutrino isocurvature velocity perturbation. 

The CMB is a powerful way to study isocurvature perturbations because 
(primordial) adiabatic and isocurvature perturbations produce very distinc- 
tive features in the CMB anisotropies ^2]- Whereas an adiabatic initial 
perturbation generates a cosine oscillatory mode in the photon-baryon fluid, 
leading to an acoustic peak at i ~ 220 (for a flat universe), a pure isocurva- 
ture initial perturbation generates a sine oscillatory mode resulting in a first 
peak at i ~ 330. The unambiguous observation of the first peak at i ~ 220 
has eliminated the possibility of a dominant isocurvature perturbation. The 
recent observation by WMAP of the CMB polarization has also confirmed 
that the initial perturbation is mainly an adiabatic mode. But this does 
not exclude the presence of a subdominant isocurvature contribution, which 
could be detected in future high-precision experiments such as Planck. 

4.7 Super-Hubble evolution 

In the case of adiabatic perturbations, there is only one (scalar) dynamical 
degree of freedom. One can thus choose either an energy density perturbation 



30 



or a metric perturbation, to study the dynamics, the other quantities being 
determined by the constraints. 

If one considers the metric perturbation \& = $ (assuming 7r = 0), one 
can combine Einstein's equations (jllOj) . with 5p = c 2 5p (for adiabatic per- 
turbations), and (jl()7j) to obtain a second-order differential equation in terms 
of \1/ only: 

+ mil + c 2 s )¥ + [2H 1 + (1 + 3c 2 s )n] ^ - k 2 c 2 ^ = 0. (133) 

Using the background Friedmann equations, the sound speed can be reex- 
pressed in terms of the scale factor and its derivatives. For scales larger 
than the sonic horizon, i.e. such that kc s <C 7i, the above equation can be 
integrated explicitly and yields 



T n 
\p = — 

a 2 



, , a 2 (H'-H 2 ) a 
a I d ^ H 2 + P 



(134) 



where a and (3 are two integration constants. 

For a scale factor evolving like a oct p , one gets 

# = —+[3pt- p -\ (135) 

p+ 1 

They are two modes: a constant mode and a decaying mode. Note that, in 
the previous subsection on the initial conditions, we eliminated the decaying 
mode to avoid the divergence when going backwards in time. 

In a transition between two cosmological phases characterized respectively 
by the scale factors a oc t pi and a oc t P2 , one can easily finds the relation be- 
tween the asymptotic behaviours of \1/ (i.e. after the decaying mode becomes 
negligible) by using the constancy of a. This gives 

* 2 = ^r*i- (136) 
Pi + 1 

This is valid only asymptotically. In the case of a sharp transition, \1/ must 
be continuous at the transition and the above relation will apply only after 
some relaxation time. For a transition radiation/non-relativistic matter, one 
finds 

9 

^mat = jQ^rad- (137) 

In practice and for more general cases, it turns out that it is much more 
convenient to follow the evolution of cosmological perturbations by resorting 
to quantities that are conserved on super-Hubble scales. A familiar example 
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of such a quantity is the curvature perturbation on uniform- energy- density 
hyper surf aces, which can be expressed in any gauge as 

( = C-H-?-. (138) 
P' 

This is a gauge-invariant quantity by definition. The conservation equation 
()113j) can then be rewritten as 

C = — ^P nad - \v 2 {E' + v), (139) 
p + r 3 

where 5P na d is the non-adiabatic part of the pressure perturbation, defined 
by 

5P nad = 5P- c 2 s 5p. (140) 

The expression ()139|) shows that ( is conserved on super-Hubble scales in the 
case of adiabatic perturbations. 

Another convenient quantity, which is sometimes used in the literature 
instead of (, is the curvature perturbation on comoving hyper surf aces, which 
can be written in any gauge as 

-n = c + -^ q . (i4i) 

It is easy to relate the two quantities ( and 1Z. Substituting e.g. 5p = 
5p c + 37iq, which follows from the definition (see (jllljl ) of the comoving 
energy density perturbation, into ()138j) . one finds 

C = -n + ^. (142) 

Using Einstein's equations, in particular (jlllj) . this can be rewritten as 

c^-^ty 2 *- <i43) 

The latter expression shows that ( and TZ coincide in the super-Hubble limit 
k < aH. 

The quantity TZ can also be expressed in terms of the two Bardeen po- 
tentials $ and Using the momentum constraint (jl08j) and the Friedmann 
equations, one finds 

TJ 

n = y-jj(y + #$) . (144) 
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In a cosmological phase dominated by a fluid with no anisotropic stress, so 
that $ = and with an equation of state P = wp with w constant, we have 
already seen that ^ is constant with time. Since the scale factor evolves like 
a oc t p with p = 2/3(1 + w), the relation (j!44j) between TZ and \l/ reduces to 

TZ = 4±^U. (145) 
3(l + u;) V ; 

In the radiation era, TZ = (3/2)\l/, whereas in the matter era, TZ = (5/3)^, 
and since TZ is conserved, one recovers the conclusion given in Eq. ()137j) . 
During inflation, w ~ — 1 and 

6 2 2 H , . 

■" + 1 = 7"-5/P (146) 

so that 

For a scalar field, the perturbed equation of motion reads 

<5ct + 3#<$<t + ( ^ + 7" J 5a = & ($ + 3*) - 2V'$. (148) 



5 Quantum fluctuations and "birth" of cos- 
mological perturbations 

In the previous section, we have discussed the classical evolution of the cosmo- 
logical perturbations. In the classical context, the initial conditions, defined 
deep in the radiation era, are a priori arbitrary. What is remarkable with 
inflation is that the accelerated expansion can convert initial vacuum quan- 
tum fluctuations into "macroscopic" cosmological perturbations (see ^3] for 
the seminal works). In this sense, inflation provides us with "natural" initial 
conditions, which turn out to be the initial conditions that agree with the 
present observations. 



5.1 Massless scalar field in de Sitter 

As a warming-up, it is instructive to discuss the case of a massless scalar 
field in a so-called de Sitter universe, or a FLRW spacetime with exponential 
expansion, a oc exp(Ht). In conformal time, the scale factor is given by 




(149) 
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The conformal time is here negative (so that the scale factor is positive) and 
goes from — oo to 0. The action for a massless scalar field in this geometry 
is given by 



S 



d x\ 



-0 



dr] d 3 x a 4 



2^ 



1-2 

2^ 



(150) 



where we have substituted in the action the cosmological metric (fTTjl . Note 
that, whereas we still allow for spatial variations of the scalar field, i.e. inho- 
mogeneities, we will assume here, somewhat inconsistently, that the geometry 
is completely fixed as homogeneous. We will deal later with the question of 
the metric perturbations. 

It is possible to write the above action with a canonical kinetic term via 
the change of variable 

u = acfi. (151) 
After an integration by parts, the action ()150|) can be rewritten as 



S 



dr] d 3 x 



u' -V« + — u 2 
a 



(152) 



The first two terms are familiar since they are the same as in the action 
for a free massless scalar field in Minkowski spacetime. The fact that our 
scalar field here lives in de Sitter spacetime rather than Minkowski has been 
reexpressed as a time- dependent effective mass 



m 



eff 



a 
a 



(153) 



Our next step will be to quantize the scalar field u by using the standard 
procedure of quantum field theory. One first turns u into a quantum field 
denoted u, which we expand in Fourier space as 



u(rj, x) 



(2tt) 3 / 2 



d 3 k {aMv)^ 3 + ^u*Me- i%£ } , (154) 



where the a) and a are creation and annihilation operators, which satisfy the 
usual commutation rules 



a fc> a k' 



4> 4' 



0. 



%> 4 



5{k-k'). 



(155) 



The function Uk{rj) is a complex time-dependent function that must satisfy 
the classical equation of motion in Fourier space, namely 



u. 



a 
a 



u k = 0, 



(156) 
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which is simply the equation of motion for an oscillator with a time- dependent 
mass. In the case of a massless scalar field in Minkowski spacetime, this 
effective mass is zero (a" /a = 0) and one usually takes Uk = {h/2k) l ^ 2 e~ tkri 
(the choice for the normalization factor will be clear below). In the case of 
de Sitter, one can solve explicitly the above equation with a" /a = 2/rf and 
the general solution is given by 

'""•''('-iJ+^ + i)- (157) 

Canonical quantization consists in imposing the following commutation 
rules on the rj =constant hypersurfaces: 

[u(r],x),u(r],x)] = [7^(77, x), Tt u (r], x 1 )] = (158) 

and 

[u(rj, x), 7f u (r], x')] = ih6(x — x'), (159) 

where tt u = 5S/5u' is the conjugate momentum of u. In the present case, 
tt u = u' since the kinetic term is canonical. 

Subtituting the expansion ()154j) in the commutator ()159|) . and using the 
commutation rules for the creation and annihilation operators ()155|) . one 
obtains the relation 

u k u' k * — u* k u' k = ih, (160) 

which determines the normalization of the Wronskian. 

The choice of a specific function Uk{ff) corresponds to a particular pre- 
scription for the physical vacuum |0), defined by 

a s |o) = 0. (i6i) 

A different choice for Uk{rf) is associated to a different decomposition into 
creation and annihilitation modes and thus to a different vacuum. 

Let us now note that the wavelength associated with a given mode k can 
always be found within the Hubble radius provided one goes sufficiently far 
backwards in time, since the comoving Hubble radius is shrinking during in- 
flation. In other words, for sufficiently big, one gets k\r]\ 3> 1. Moreover, 
for a wavelength smaller than the Hubble radius, one can neglect the influ- 
ence of the curvature of spacetime and the mode behaves as in a Minkowski 
spacetime, as can also be checked explicitly with the equation of motion ()156j) 
(the effective mass is negligible for A;|?7| ^> 1). Therefore, the most natural 
physical prescription is to take the particular solution that corresponds to 
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the usual Minkowski vacuum, i.e. Uk ~ exp(— ikr)), in the limit k\rj\ ^ 1. In 
view of ()157|) . this corresponds to the choice 




(162) 



where the coefficient has been determined by the normalisation condition 
(jlfiOj) . This choice, in the jargon of quantum field theory on curved space- 
times, corresponds to the Bunch-Davies vacuum. 

Finally, one can compute the correlation function for the scalar field 4> hi 
the vacuum state defined above. When Fourier transformed, the correlation 
function defines the power spectrum V^ik\. 

(010(^)0(^)10) = J d 3 k e **(*-a)^M. (163) 

Note that the homogeneity and isotropy of the quantum field is used implic- 
itly in the definition of the power spectrum, which is "diagonal" in Fourier 
space (homogeneity) and depends only on the norm of k (isotropy). In our 
case, we find 

2n 2 k~ 3 V, = (164) 

which gives in the limit when the wavelength is larger than the Hubble radius, 
i.e. k\r]\ <C 1, 

V^k) ~ % (j^j (k < aH) (165) 

Note that, in the opposite limit, i.e. for wavelengths smaller than the Hubble 
radius (k\rj\ ^> 1), one recovers the usual result for fluctuations in Minkowski 
vacuum, V^{k) = Ti{k/2na) 2 . 

We have used a quantum description of the scalar field. But the cos- 
mological perturbations are usually described by classical random fields. 
Roughly speaking, the transition between the quantum and classical (al- 
though stochastic) descriptions makes sense when the perturbations exit the 
Hubble radius. Indeed each of the terms in the Wronskian ()160)) is roughly 
of the order h/2(kr]) 3 in the super-Hubble limit and the non-commutativity 
can then be neglected. In this sense, one can see the exit outside the Hubble 
radius as a quantum-classical transition, although much refinement would be 
needed to make this statement more precise. 

5.2 Quantum fluctuations with metric perturbations 

Let us now move to the more realistic case of a perturbed inflaton field living 
in a perturbed cosmological geometry. The situation is more complicated 
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than in the previous problem, because Einstein's equations imply that scalar 
field fluctuations must necessarily coexist with metric fluctuations. A correct 
treatment, either classical or quantum, must thus involve both the scalar field 
perturbations and metric perturbations. 

In order to quantize this coupled system, the easiest procedure consists in 
identifying the true degrees of freedom, the other variables being then derived 
from them via constraint equations. As we saw in the classical analysis of 
cosmological perturbations, there exists only one scalar degree of freedom in 
the case of a single scalar field, which we must now identify. 

The starting point is the action of the coupled system scalar field plus 
gravity expanded up to second order in the linear perturbations. Formally 
this can be written as 

S[0+50, g, v = g^+h^] = [0, g» v ]+S^ [50, h^, 0, g^+S® [50, h 

(166) 

where the first term contains only the homogeneous part, contains 
all terms linear in the perturbations (with coefficients depending on the ho- 
mogeneous variables), and finally contains the terms quadratic in the 
linear perturbations. When one substitutes the FLRW equations of motion 
in (after integration by parts), one finds that vanishes, which is 
not very surprising since this is how one gets the homogeneous equations of 
motion, via the Euler-Lagrange equations, from the variation of the action. 

The term is the piece we are interested in: the corresponding Euler- 
Lagrange equations give the equations of motion for the linear perturbations, 
which we have already obtained; but more importantly, this term enables us 
to quantize the linear perturbations and to find the correct normalization. 

If one restricts oneself to the scalar sector, the quadratic part of the action 
depends on the four metric perturbations A, B, C, E, as well as on the scalar 
field perturbation 50. After some cumbersome manipulations, by using the 
FLRW equations of motion, one can show that the second-order action for 
scalar perturbations can be rewritten in terms of a single variable 



v = a 50 -^-C ), (167) 




which is a linear combination mixing scalar field and metric perturbations. 
The variable v represents the true dynamical degree of freedom of the system, 
and one can check immediately that it is indeed a gauge-invariant variable. 

In fact, v is proportional to the comoving curvature perturbation defined 
in ()141j) and which, in the case of a single scalar field, takes the form 

TZ = -C + — 50. (168) 
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Note also that, if one can check a posteriori that v is the variable describ- 
ing the true degree of freedom by expressing the action in terms of v only 
(modulo, of course, a multiplicative factor depending only on homogeneous 
quantities: the v defined here is such that it gives a canonical kinetic term 
in the action), one can identify v in a systematic way by resorting to Hamil- 
tonian techniques, in particular the Hamilton- Jacobi equation |15j . 

With the variable v, the quadratic action takes the extremely simple form 



S v = - J drj d 3 x 



v' 2 + diV&V + — v 2 

z 



(169) 



with 



z = o^. (170) 

This action is thus analogous to that of a scalar field in Minkowski spacetime 
with a time-dependent mass. One is thus back in a situation similar to 
the previous subsection, with the notable difference that the effective time- 
dependent mass is now z"/z, instead of a" /a. 

The quantity we will be eventually interested in is the comoving curvature 
perturbation TZ, which is related to the canonical variable v by the relation 

v = -zK. (171) 

Since, by analogy with ()164|) . the power spectrum for v is given by 

2ir 2 k~ 3 V v (k) = \v k \ 2 , (172) 

the corresponding power spectrum for 1Z is found to be 

2n 2 k- 3 V n (k) = (173) 
z z 

In the case of an inflationary phase in the slow-roll approximation, the 
evolution of <fi and of H is much slower than that of the scale factor a. 
Consequently, one gets approximately 

z" a!' 

— ~ — , (slow - roll) (174) 

z a 

and all results of the previous section obtained for u apply directly to our 
variable v in the slow-roll approximation. This implies that the properly 
normalized function corresponding to the Bunch-Davies vacuum is approxi- 
mately given by 
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In the super-Hubble limit fc|77| <C 1 the function behaves like 

hi h aH 

Wfc * "V ^ ^'V 2* X' (176) 

where we have used a ~ —\/{Hrj). Consequently, on scales larger than the 
Hubble radius, the power spectrum for 71 is found, combining l| 173)1 . (J17()j) 
and 1)175)1 . to be given by 

where we have reintroduced the cosmic time instead of the conformal time. 
This is the famous result for the spectrum of scalar cosmological perturba- 
tions generated from vacuum fluctuations during a slow-roll inflation phase. 
Note that during slow-roll inflation, the Hubble parameter and the scalar 
field velocity slowly evolve: for a given scale, the above amplitude of the per- 
turbations is determined by the value of H and <j) when the scale exited the 
Hubble radius. Because of this effect, the obtained spectrum is not strictly 
scale-invariant. 

It is also instructive to recover the above result by a more intuitive deriva- 
tion. One can think of the metric perturbations in the radiation era as result- 
ing from the time difference for the end of inflation at different spatial points 
(separated by distances larger than the Hubble radius), the shift for the end 
of inflation being a consequence of the scalar field fluctuations 5<p ~ H/2tt. 
Indeed, 

^ rad ~ — ~ H5t, (178) 



and the time shift is related to the scalar field fluctuations by St ^ 
which implies 

H 2 

^rad ~ -j", (179) 

which agrees with the above result since during the radiation era 1Z — (3/2) ^ 
(see Eq. (1145)) ). It is also worth noticing that, during inflation, in the case of 
the slow-roll approximation, the term involving C in the linear combination 
(Jl(i7j) defining v is negligible with respect to the term involving 6(f). One 
can therefore "ignore" the role of the metric perturbations during inflation 
in the computation of the quantum fluctuations and consider only the scalar 
field perturbations. But this simplification is valid only in the context of 
slow-roll approximation. It is not valid in the general case, as can be verified 
for inflation with a power-law scale factor. 
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5.3 Gravitational waves 



We have focused so far our attention on scalar perturbations, which are the 
most important in cosmology. Tensor perturbations, or primordial gravita- 
tional waves, if ever detected in the future, would be a remarkable probe 
of the early universe. In the inflationary scenario, like scalar perturbations, 
primordial gravitational waves are generated from vacuum quantum fluctu- 
ations. Let us now explain briefly this mechanism. 

The action expanded at second order in the perturbations contains a 
tensor part, which given by 

Sf = / drj d 3 x a 2 ^d,E)d u El (180) 

where 7]^ denotes the Minkoswki metric. Apart from the tensorial nature of 
Ej, this action is quite similar to that of a scalar field in a FLRW universe 
(|150p . up to a renormalization factor l/yWhiG. The decomposition 

E I if^^M^k-A)^ (181) 

where the e* (k; A) are the polarization tensors, shows that the gravitational 
waves are essentially equivalent to two massless scalar fields (for each polar- 
ization) <p\ = m P E\/2. 

The total power spectrum is thus immediately deduced from ()164|) : 

p » = 2x i xf! (s) 2 - < 182 > 

where the first factor comes from the two polarizations, the second from the 
renormalization with respect to a canonical scalar field, the last term being 
the power spectrum for a scalar field derived earlier. In summary, the tensor 
power spectrum is 

2h ( H \ 2 . . 

V 9 = — \ — ) > 183 

vr 2 \m P J k=aH 

where the label recalls that the Hubble parameter, which can be slowly evolv- 
ing during inflation, must be evaluated when the relevant scale exited the 
Hubble radius during inflation. 

5.4 Power spectra 

Let us rewrite the scalar and tensor power spectra, respectively given in (j 177)1 
and l| 183)1 . in terms of the scalar field potential only. This can be done by 
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using the slow-roll equations (j49H50j) . One finds for the scalar spectrum 



with subscript meaning that the term on the right hand side must be evalu- 
ated at Hubble radius exit for the scale of interest. The scalar spectrum can 
also be written in terms of the first slow- roll parameter defined in ()52|). in 
which case it reads 



1 



24tt 2 \m%e v J k=aH 
From the observations of the CMB fluctuations, 

r]l 2 = — U ( ~~T7^\ - 5 x 10 ~ 5 - ( 186 ) 

If ey is order 1, as in chaotic models, one can evaluate the typical energy 
scale during inflation as 

V l/i ~ 10~ 3 m P ~ 10 15 GeV. (187) 
The tensor power spectrum, in terms of the scalar field potential, is given 

by 

P - = ^S) M - (188) 

The ratio of the tensor and scalar amplitudes is proportional to the slow-roll 
parameter ey: 

r = ^ = 16e v . (189) 

r-R, 

The scalar and tensor spectra are almost scale invariant but not quite 
since the scalar field evolves slowly during the inflationary phase. In order to 
evaluate quantitatively this variation, it is convenient to introduce a scalar 
spectral index as well as a tensor one, defined respectively by 

ns{k) - 1 = -^r 1 ' " tW = ^w • (190) 

One can express the spectral indices in terms of the slow-roll parameters. 
For this purpose, let us note that, in the slow-roll approximation, dink = 
dln(aH) ~ din a, so that 

* * V ' ~-mlt. (191) 



dlna H 3H 2 P V 
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where the slow-roll equations (j49H5U|) have been used. Therefore, one gets 



n s (k) -1 = 2r] V -6e v , (192) 

where ey and rjv are the two slow-roll parameters defined in ()52j) and t)54j). 
Similarly, one finds for the tensor spectral index 

n T (k) = -ley. (193) 

Comparing with Eq. ()189|) . this yields the relation 

r = -8n T , (194) 

the so-called consistency relation which relates purely observable quantities. 
This means that if one was able to observe the primordial gravitational waves 
and measure the amplitude and spectral index of their spectrum, a rather 
formidable task, then one would be able to test directly the paradigm of 
single field slow-roll inflation. 

Finally, let us mention the possibility to get information on inflation 
from the measurement of the running of the spectral index. Introducing 
the second-order slow-roll parameter 

A V'V" 

iv = m%-^, (195) 

the running is given by 

dn s „ . 2 



dink 



-24e 2 v + We v r] V - 2£ v . (196) 



As one can see, the amplitude of the variation depends on the slow-roll 
parameters and thus on the models of inflation. 



5.5 Conclusions 

To conclude, let us mention the existence of more sophisticated models of 
inflation, such as models where several scalar fields contribute to inflation. 
In contrast with the single inflaton case, which can generate only adiabatic 
primordial fluctuations, because all types of matter are decay product of 
the same inflaton, multi-inflaton models can generate both adiabatic and 
isocurvature perturbations, which can even be correlated [TH] . 

Another recent direction of research is the possibility to disconnect the 
fluctuations of the inflaton, the field that drives inflation, from the observed 
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cosmological perturbations, which could have been generated from the quan- 
tum fluctuations of another scalar field [T7j . 

From the theoretical point of view an important challenge remains to 
identify viable and natural candidates for the inflaton field (s) in the frame- 
work of high energy physics, with the hope that future observations of the 
cosmological perturbations will be precise enough to discriminate between 
various candidates and thus give us a clue about which physics really drove 
inflation. 

Alternative scenarios to inflation can also been envisaged, as long as 
they can predict unambiguously primordial fluctuations compatible with the 
present observations. In this respect, one must emphasize that the cosmolog- 
ical perturbations represent today essentially the only observational window 
that gives access to the very high energy physics, hence the importance for 
any early universe model to be able to give firm predictions for the primordial 
fluctuations it generates. 
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